Вестник Донского государственного технического университета 


2016, №4(87), 118-125 





ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ 
ТЕХНИКА И УПРАВЛЕНИЕ 
ПМЕОВМАТТОМ ТЕСНМОГОСУ, СОМРОТЕКВ 
ЭСТЕМСЕ, АХО МАМАСЕМЕМТ 


УДК 539.182; 519.632 





РОГ 10.12737/22156 


Применение численных фундаментальных решений в методе точечных источников поля" 


С. Ю. Князев', Е. Е. Щербакова?” 


Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 


АррИсайоп оЁ Фе питечсаПу оматед Рапдашениа! зобопз и Ше йе ройи-зоигсе шефо@””” 


5. Уи. Кпуахеу', Е. Е. ЗвспегфаКота?”” 


1.2 роп {ще Тесыиса1 Ошщуегзиу, Возюу-оп-Ооп, Киззап Еедеганоп 


Целью работы является получение интегрального уравнения, 
с помощью которого, используя известное фундаментальное 
решение другого уравнения, возможно численным методом 
найти фундаментальное решение линейного уравнения эллип- 
тического типа. Вводится понятие численного фундаменталь- 
ного решения (ЧФР). Полученные таким образом численные 
фундаментальные решения (ЧФР) могут быть использованы 
при решении краевых задач для уравнений эллиптического 
типа различной размерности с помощью метода точечных 
источников поля (МТИ). Результатом работы является созда- 
ние эффективного численного метода решения краевых задач 
с использованием ЧФР. Это позволяет расширить круг реша- 
емых с помощью МТИ задач. Таким образом, МТИ выступает 
в качестве универсального численного метода при решении 
краевых задач для линейных уравнений эллиптического типа. 
Особенно эффективно применение предложенного способа 
при решении трехмерных задач Дирихле для уравнений со 
сферически симметричными фундаментальными решениями. 
В качестве тестовой задачи предложенным способом решено 
уравнение Шредингера для одномерного квантового осцилля- 
тора. Показано, что, используя фундаментальные решения 
уравнения Шредингера, полученные численно, удается найти 
собственные значения и собственные функции квантового 
осциллятора. Найденные собственные функции осциллятора 
соответствуют известным аналитическим решениям кванто- 
вой задачи. В качестве другого тестового примера решается 
двумерная краевая задача для уравнения Гельмгольца. В этом 
случае предварительно находится численное фундаменталь- 
ное решение для уравнения Гельмгольца. Вычислены зависи- 
мости погрешности численного решения от числа узлов в 
области решения задачи. На основании полученных результа- 
тов делается вывод о перспективности предложенного чис- 
ленного метода. 


Ключевые слова: фундаментальное решение, метод фунда- 
ментальных решений, метод точечных источников, метод 
интегрированных источников, дискретные источники. 
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Введение. Метод точечных источников поля (МТИ) является одним из эффективных методов моделирования 
физических полей (например, электрических и магнитных) в технических (в том числе электромеханических) устрой- 
ствах [1-4]. Для указанного метода характерна высокая точность численного решения и чрезвычайная простота ком- 
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пьютерной реализации [5-9]. Наилучшие результаты получены при использовании МТИ для моделирования физиче- 
ских полей, описываемых однородными линейными уравнениями эллиптического типа с известными фундаменталь- 
ными решениями, задаваемыми аналитически. Такими уравнениями являются уравнение Лапласа, уравнение Гельм- 
гольца, бигармоническое уравнения, некоторые другие типы уравнений. Применение МТИ в этом случае позволило 
решить значительное число прикладных задач по моделированию, например: стационарных электрических, магнит- 
ных [8-10], тепловых, концентрационных полей [11-14], полей упругих напряжений [15—17]. МТИ успешно применя- 
ется также при численном решении краевых задач для неоднородных уравнений, таких как уравнение Пуассона [18— 
19], неоднородное уравнение Гельмгольца [19—21]. Однако во всех случаях использования МТИ предполагается из- 
вестным фундаментальное решение соответствующего уравнения математической физики. Это резко ограничивает 
круг решаемых с помощью МТИ задач. Тем не менее, и в этом случае возможно решение краевых задач с помощью 
МТИ, если предварительно найти численные значения фундаментальных решений при определенных значениях па- 
раметров. Назовем фундаментальное решение, заданное численно, численным фундаментальным решением (ЧФР). 
Ниже описан метод нахождения ЧФР для линейных уравнений эллиптического типа и показывается возможность ис- 
пользования этих решений в МТИ. 

Интегральное уравнение для дополнительной составляющей фундаментального решения. Пусть фун- 


даментальное решение $(г) линейного уравнения эллинтического типа [0 =0 известно. Требуется найти фунда- 
ментальное решение С(г) уравнения типа (+10 =0. Здесь Г и { — линейные операторы эллиптического типа. 
Представим искомое фундаментальное решение С (г) в виде 
(к) = в(г)и(). (р 
Здесь и(г) — искомая функция. Назовем &(г) основной составляющей, а и(г) — дополнительной составляющей 
фундаментального решения С(х). 
В соответствии с определением фундаментального решения справедливо соотношение 
(2+06=5(). 
Подставим в это выражение (1). В результате получим 
(2+1) (8 +и) = 5(,) : 
или [8+1 -1е+ш =56(х). (2) 
Так как Ге = 5(,) (2 (г) есть фундаментальное решение уравнения Г. =0 ), то уравнение (2) запишем в виде 
Ги =—ю-—Ш. (3) 
Правую часть этого уравнения будем рассматривать как плотность заряда р(г) неоднородного уравнения 


Ги =р (г) . Тогда частное решение неоднородного уравнения (3) можно записать в виде 


и(г) = [8 (к, В)р(®)ао» 
или в виде 


ва В АВЕ . (4) 


Областью интегрирования в (4) может быть шар радиусом Ко, который заключает в себя область @® задачи, 
решаемой с помощью искомого фундаментального решения. 

Уравнение (4) допускает численное решение. В результате находятся приближенные значения для ЧФР в уз- 
ловых точках области ©) > ®. Эти значения можно использовать при решении задачи Дирихле с помощью, например, 
МТИ. При решении задачи Неймана или краевой задачи третьего рода наряду со значениями фундаментального ре- 
шения необходимо иметь возможность вычисления нормальных производных фундаментального решения на границе 
области © или частных производных по координатам. Для вычисления частных производных можно воспользоваться 
соотношением (4). Например, продифференцировав (4) по координате х, получим 

ди (г) де (г, К) 
дх 1 дх 





(и(®)+3(К))аОьв . (5) 


Аналогичное соотношение можно получить для частной производной дополнительной составляющей фунда- 
ментального решения по координате у и по координате < (при решении трехмерных задач). Используя найденные зна- 


чения функции и(г) в узлах области © с помощью численного интегрирования в правой части соотношения (5) мож- 


НО ВЫЧИСЛИТЬ приближенные значения частной производной дополнительной составляющей фундаментального реше- 
ния в любой точке области © и использовать эти значения при численном решении краевой задачи. 
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Фундаментальные решения уравнения Шредингера. В качестве простейшего примера найдем фундамен- 
тальные решения одномерного уравнения Шредингера, описывающего колебания квантового осциллятора. Пусть ча- 


стица массой т совершает упругие колебания с частотой © и энергией Е. Уравнение Шредингера запишем в виде [22] 





2 22 
2 
ы У Ву =“ > ху, 
ах й 
| то 
а после замены переменных ит. —›х запишем как 
2 
Ее = ху ‚ или как 
4? То 
2 
\ 2 
чер = ху, (6) 
ах 
2 
гдее=— ЕЁ. 
йо 
2 
Здесь операторами Г и [, фигурирующими в (2)-(4), будем считать Ё = и = —х?. Как известно, фунда- 


Ах 


4? и (№ |) 
ментальное решение уравнения о еу=0 равно 8 (х) = 


/х 2 


[23]. Классически допустимая область движе- 


ния осциллятора ограничивается отрезком [-\е, \е] › ПОЭТОМУ В ее области решения для квантового осциллято- 


ра можно взять отрезок = Уее |=, И, а областью интегрирования в (4) считать, например, отрезок 


О = Аа, 1,5/е] =[-а,4]. Тогда уравнение (4) запишется как 


о (Меж) х?и(х)ах =. р (Мех (| х?ах | (т) 


2\/е Де 


Обозначим интеграл в правой части интегрального уравнения (7) как /(х). В данном случае функция }(х) вы- 
числяется аналитически. В более общем случае, если исходное дифференциальное уравнение имеет другой вид, функ- 
цию / (>) находят численно. 

Для численного решения уравнения (7) интеграл в левой части соотношения (7) представляют в виде инте- 
гральной суммы, используя тот или иной метод численного интегрирования. Если, например, используется простей- 
ший метод численного интегрирования, метод прямоугольников, то отрезок © разбивается на № элементарных отрез- 


КОВ. 
Обозначим значения искомой функции и(х) в некоторых, например, в средних точках х; элементарных от- 


резков как и;, где 1— номер элементарного отрезка. В результате для величин иИ,, будет получена система линейных 


алгебраических уравнений вида 





1, =} 
(И: =. Ь. = ‚а, = Й. , р 
Хахи, ‚› ГДЕ 0; 1(х;) ау т 


Здесь й; — длина элементарного отрезка с номером /. 
После того как величины и; найдены, значения искомого фундаментального решения в узлах х; , т. е. ЧФР вы- 


числяются с помощью формулы (1): 
С(х)=8(х,)+щ,. 
Для дальнейшего удобно (но не обязательно), используя величины С (х; )‚ получить интерполяционную фор- 
мулу для фундаментального решения С (х). Теперь можно приступить к нахождению решения уравнения (5) с помо- 
щью МТИ. Для этого справа и слева от области решения задачи, т. е. отрезка ®, в точках с координатами 5, =-[-бх и 


<, =[+0х размещаются заряды 4; и 4, , моделирующие искомое поле \\(х) в области ®. Здесь бх — удаленность за- 
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рядов, моделирующих поле \у(х), от границ области в. Разумеется, должно выполняться условие [-+бх < 4 ‚т. е. точки 
< и <, не должны выйти за пределы области ©. Искомое решение представляем как суперпозицию полей двух то- 
чечных зарядов, т. е. в виде соотношения 
у(х) = 4,6 (1+х+5х)+а,С(1-х+5х). (8) 
Для окончательного решения задачи необходимо найти величины зарядов 4; и 4, ‚ используя два каких-либо 
дополнительных условия, накладываемых на волновую функцию \у(х). Очевидно, что физический смысл имеют толь- 
ко четные и нечетные волновые функции. Отсюда следует, что заряды 4, и 4, либо равны друг другу (а, =4,=а), 
либо имеют противоположные знаки (4, =-4, =4). Величину заряда 4 следует выбрать так, чтобы удовлетворить 


условию нормировки ВОЛНОВОЙ функции, например, потребовав, чтобы максимальное по абсолютной величине значе- 
ние волновой функции в области решения < равнялось единице. 


Известно [22], что решение уравнения (6) имеет реальный физический смысл не при любых значениях энер- 





йо 1 
гии А, т. е. параметра е, а только при Ё = ь е=й®| п+ ой где квантовое число и может принимать только целочис- 


ленные значения, начиная с нуля. Только при этих значениях квантового числа волновая функция на бесконечности 
обращается В НОЛЬ. Очевидно, что при этом на границах области решения волновая функция будет иметь значения 


заметно меньшие, чем при других, не целочисленных значениях числа п. 


Назовем величину К =1/ (№(-] + (1) критерием задачи. На рис. 1 представлен график зависимости кри- 


терия задачи Кот энергии частицы, точнее, от числа п, которое изменяется в пределах от 0 до 10. Находились четные 


решения уравнения Шредингера. Вычисления проводились при числе узлов в области © М = 100. 


Параметр К 





0 20 40 60 80 100 
10 


Рис. 1. Зависимость параметра К =1/ (м (-1) + № (1 )) от квантового числа и 


Из рис. 1 видно, что максимумы на кривой соответствуют значениям квантового числа п, близким к 0, 2, 4, 6, 
8 и 10. Это в точности согласуется с известными результатами квантовой теории для собственных значений энергии и 
для состояний осциллятора с четными волновыми функциями. Соответствующая зависимость К = К(п), полученная 
для нечетных волновых функций, имеет максимумы при нечетных целочисленных значениях квантового числа и. 
Решение уравнения Шредингера (5) в аналитическом виде [22]: 
2 


х 
ци, (х) = Аехр - Н, (Хх): (9) 
где Н,„ (х) — полином Эрмита степени п, сравнивалось с волновой функцией \у(х), рассчитанной по формуле (8), и 


погрешность численного метода оценивалась с помощью формулы 


== пах |, (х)- (х) . 


Вестник Донского государственного технического университета 2016, №4(87), 118-125 





При вычислении погрешности = волновая функция (9) нормировалась так же, как и волновая функция \у(х), 


т. е. нормировочный множитель А в (9) выбирался таким, чтобы максимальное по абсолютной величине значение вол- 
Новой функции я (х) в области решения о равнялось единице. Очевидно, что точность численного решения зависит 


от числа узлов в области О, т. е. от числа №. На рис. 2 представлены такие зависимости, полученные при значениях 


кантового числа и = 0 (жирная сплошная линия); и = 2 (тонкая сплошная линия); п = 4 (пунктирная линия). 
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Рис. 2. Зависимость погрешности численного решения для ВОЛНОВОЙ функции от числа узлов М№в области решения для разных зна- 
чений квантового числа и: и = 0 (жирная сплошная линия); п = 2 (тонкая сплошная линия); и = 4 (пунктирная линия) 

Таким образом, применение численного фундаментального решения позволяет не только найти волновые 
функции, т. е. решить уравнение Шредингера, но и установить собственные значения энергии. Точность полученных 
результатов легко регулируется путем изменения числа узлов М в области решения задачи. Снизить погрешность ре- 
шения можно также, используя более точный метод численного интегрирования. 

Тестовый пример решения двумерной краевой задачи. В качестве второй тестовой задачи решалась задача 
Дирихле для уравнения Гельмгольца 

ДИ +^И =0 
в круговой области ® радиусом го. При использовании МТИ заряды, моделирующие искомое поле, располагались на 
вспомогательной окружности радиусом К. = К-го, где постоянная К > 1 определяет удаленность моделирующих заря- 


дов от границ области ®. При решении данной тестовой задачи применялись ЧФР для уравнения Гельмгольца, вычис- 


в 1 
ленные с использованием фундаментальных решений для уравнения Лапласа з (г, К ) = и. |" _—Ю | . Таким образом, в 
п 


данной задаче [=А и [=^. В качестве области интегрирования © в (4) использовался круг радиусом Ко. На гео- 


метрические параметры задачи накладывалось условие /) А, < К. Моделировалось поле Гельмгольца точечного 


заряда, помещенного в точку р = (хо,0) на оси Х. Известные потенциалы этого поля сравнивались с вычисленными с 
помощью МТИ при использовании ЧФР значениями, и вычислялась относительная среднеквадратичная погрешность 
& решения. При численном интегрировании в уравнении (4) использовались значения подинтегральных функций в 
узлах равномерной прямоугольной сетки. Количество узлов М сетки, находящихся внутри или на границе области ©, 
является одним из важнейших параметров задачи, определяющим точность ее решения. Часть этих узлов в количестве 
№, находится внутри области решения задачи ®. Другим важным параметром является количество моделирующих 
зарядов №., равномерно располагаемых на вспомогательной окружности. На рис. 3 представлены зависимости средне- 
квадратичной погрешности = от количества узлов №, в области ®. Вычисления производились при следующих значе- 


ниях параметров: №, = 50; го =2; К = 1,3; К. = Кто; Ко =5. 
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Рис. 3. Зависимость среднеквадратичной погрешности от количества узлов №, в области решения задачи © для разных значений 
параметра ^: ^. = -0,5 (черная пунктирная линия); /. = —1,0 (черная сплошная линия); /, = 0,5 (красная пунктирная линия); 
\. = 1,0 (красная сплошная линия) 


Кривые, представленные на рис. 3 черным цветом, получены при отрицательных значениях параметра /: 
пунктирной линии соответствует ^. = -0,5; сплошной линии соответствует /. = —1,0. Кривые, представленные на рис. 3 
красным цветом, получены при положительных значениях параметра /: пунктирной линии соответствует 4. = 0,5; 
сплошной линии соответствует /, = 1,0. Из данных рис. 3 видно, что при решении краевой задачи предложенным спо- 
собом легко достигается решение с относительной среднеквадратичной погрешностью порядка 10-10. Этот ре- 
зультат можно считать вполне приемлемым. Кроме того, из рис. 3 видно, что, как и следовало ожидать, погрешность 
результата убывает с ростом количества узлов №, в области ©. Можно ожидать, что при использовании более точного, 
по сравнению с методом прямоугольников, метода численного интегрирования в выражении (4) можно получить бо- 
лее точное численное решение краевой задачи и более быстрое снижение погрешности с ростом числа узлов М,.. 


Выводы Рассмотренные выше тестовые примеры подтверждают возможность использования ЧФР при реше- 
нии краевых задач с помощью МТИ. Особенно эффективно применение предложенного способа при решении задач 
Дирихле для уравнений со сферически симметричными фундаментальными решениями. Таким свойством обладают, 
например, фундаментальные решения для уравнений Лапласа и Гельмгольца. В этом случае фундаментальное реше- 
ние зависит только от величины аргумента, т. е. от длины радиус-вектора г, а не от его ориентации. Поэтому в урав- 
нении (4) можно предварительно произвести интегрирование по угловым координатам и получить одномерное инте- 


гральное уравнение для численного нахождения функции и(г) и вычисления ЧФР, которое в этом случае представля- 


ет собой одномерный массив чисел. Так как после нахождения ЧФР вычисление производных для фундаментальных 
решений легко производится с помощью соотношений подобных (5), то это делает возможным решение задач Нейма- 
на и краевых задач третьего рода. Следует отметить также, что ЧФР, полученные для некоторой области ©, могут ис- 
пользоваться в МТИ при решении различных краевых задач для различных областей осо. 
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